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We present near ultraviolet (NUV:1750 - 2800A) and far 
ultraviolet (FUV: 1350 - 1750A) light-curves for flares on 
4 nearby dMe-type stars (GJ 3685A, CR Dra, AF Psc and 
SDSS J084425.9+513830.5)observed with the GALEX satel- 
lite. Taking advantage of the time-tagged events recorded with 
the GALEX photon counting detectors, we present high tem- 
poral resolution (< 0.01 sec) analysis of these UV flare data. 
A statistical analysis of 700 seconds of pre-fiare quiescence 
data for both CR Dra and SDSS J084425.9+513830.5 failed 
to reveal the presence of significant micro-flare activity in time 
bins of 0.2, 1 and 10 second intervals. Using an appropriate 
differential emission measure for both the quiescent and flar- 
ing state, it is possible to reproduce the observed FUV:NUV 
flux ratios. A major determinant in reproducing this flux ra- 
tio is found to be the value of plasma electron density during 
the flare. We also searched the count rate data recorded dur- 
ing each of the four flare events for periodicity associated with 



^ ■ magneto-hydrodynamic oscillations in the active region coro- 
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nal loops. Significant oscillations were detected during the flare 
events observed on all 4 stars, with periodicities found in the 30 
,— ' \ to 40 second range. Flare oscillations with this periodicity can 
,— ' be explained as acoustic waves in a coronal loop of length of 
ps 10 9 cm for an assumed plasma temperature of 5 — 20 x 10 6 K. 
This suggests a loop length for these M-dwarf flares of less than 
Tj" I 1 /10 th of the stellar radii. We believe that this is the first de- 
tection of non-solar coronal loop flare oscillations observed at 
^ \ ultraviolet wavelengths. 
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1. Introduction 

Although flare stars were detected early in the last century, it 
was not until the late 1940's that flare star research attracted 
wide-spread attention with the observation by Joy & Humason 
J1949i of a four magnitude eruptive event. Since then there 
have been many hundreds of publications concerning these ac- 
tive stellar sources. The first coincidence between optical and 
radio activity was reported by Lovell &Solomon ( 1966), with 
the first X-ray flare being observed by Heise et al. (1975 1. 
Surveys show that low mass M-dwarf (dMe) stars account for 
more than 75% of the stellar population in the solar neighbor- 
hood and, due to the strong magnetic fields that cover most of 
their stellar disks, many of them exhibit significant levels of 
coronal activity. In a recent survey of 8000 late-type dwarfs 
contained in the Sloan Digital Sky Survey (SDSS) catalog, 
West et al. (2004) have found that > 50% of stars of spec- 
tral class M4 to M9 have high levels of magnetic activity, as 
inferred by the presence of the Ha emission line in their stel- 
lar spectra. The most easily observable manifestation of this 
activity is that of stellar flare eruptions, as typified by the star 
UV Ceti, which involves random outbursts in which the short- 
term stellar brightness increases significantly on time-scales of 
seconds to hours at X-ray, UV, visible and radio wavelengths 
(Schmitt et al. [T993l Phillips et al. [T988l Stepanov et al. [T995V 
Although the detailed physics associated with the formation 
of stellar flares is still uncertain (Haisch, Strong & Redono 
I199U . it is clear that these eruptions are linked to the mag- 
netic heating processes occurring in stellar coronae. We note 
that the soft X-ray emission from flares (which is indicative of 
gas with temperatures of > 10 6 K) is dominated by coronal line- 
emission, while the UV radiation during flare events originates 
in the chromosphere and transition region and is characterized 
by emission from a gas with a temperature of ~ 10 5 K. There 
is still much debate as to the relative contributions from emis- 
sion line and/or UV continuum flux to the total heating budget 
in flare stars (Robinson et al. 120011 Hawley et al 2003 Giidel 
et al 2003 1. Recent EUV and X-ray observations suggest that 



flares statistically contribute a significant fraction of the over- 
all (X-ray) heating of the coronae on all magnetically active 
stars (Giidel et al. 120031 Arzner k Gtidel 12004k In particular, 
the recent XMM — Newton satellite observations of dMe star 
flares by Mitra-Rraev et al. ( 2005a I have noted a correlation be- 
tween the observed ultraviolet (UV) flare energy and the corre- 
sponding increase in X-ray luminosity, with the UV flares tem- 
porally preceding the peaks in X-ray flux. These observations 
would seem to support the predictions for the Neupert Effect 
(Neupert ll968> . as outlined in the chromospheric evaporation 
model for solar flares (Antonucci, Gabriel & Dennis 1984}- In 
this picture, during the impulsive phase of a stellar flare, the 
UV/optical emission due to accelerated electrons that gyrate 
downward along the magnetic field lines and impact the chro- 
mosphere, should precede the more slowly evolving X-rays that 
are emitted by the heated plasma. 

Although flare events have been routinely recorded on 
many M-dwarfs at radio( Jackson, Kundu & White 1 989j>, op- 
tical (Gunn et al. [T994t . Ultra-violet (Robinson et al. I200T1 . 
Extreme Ultra-violet (Audard et al. 2000 1 and X-ray (Marino, 
Micela & Peres 2000) wavelengths, such observations have 
generally suffered from both poor temporal-resolution and the 
inability to track a large event throughout the entire pre-flare 
to post-flare period. However, Beskin et al. (1988) recorded 
emission structure on time-scales of 0.3 - 0.8 seconds in high- 
temporal (0.3 milli-second) U-band visible observations of 
flare events on several bright UV Ceti-type stars. Using the 
scale height of the atmosphere and the velocity of the shock 
led them to the conclusion that optical flares on M-dwarfs were 
due to thermal phenomena. We note that U-band observations 
of both solar and stellar flares have been widely used as a 
proxy for non-thermal hard X-ray flare emission (Neidig[l989 
Hawley et al.[T333}. 

High (< 1 sec) temporal resolution observations of ultra- 
violet flares can provide detailed information on microflaring 
activity and also determine whether wave activity is present. 
Both of these are thought to be important in the heating of 
coronae (Robinson, Carpenter & Percival fl999l Roberts 2004). 
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Star Name R.A.(J2000) Dec.(J2000) Survey Mode MAST Identifier Exposure (sec) 

GJ 3685A 11:47:41 +00:15:12 MIS MISDR1 - 1 3062-0283-0003 1244 

AFPsc 23:31:44 -02:45:12 GI* GI1-067026-PGC7 1626-0002 1705 

CRDra 16:17:05 +55:16:09 DIS ELAISN 1-07-0004 1650 

SDSSJ084425 08:44:26 +51:38:31 MIS MISDR 1-03333-0447-0001 1669 

* pointed Guest Investigator observation 



Furthermore, the ability to record emission throughout an en- 
tire energetic flare event can provide insights into the nature of 
the pre-flare and post-flare quiescent phases and how their sta- 
tistical small-scale flux variations may be related to the onset 
of a major flare event ( Gudel et al. l20U4l The NASA GALEX 
satellite (Martin et al. I2005> . launched in 2003, has recently 
been shown to be an excellent observational platform for the 
serendipitous detection and high time-resolved observation of 
stellar flares in two ultraviolet (UV) photometric bands. For 
example, an extremely large flare event was observed on the 
dM4e star GJ 3685A on April 20th, 2004 that revealed an in- 
crease in stellar brightness of > 10 UV magnitudes in a period 
of < 200 seconds (Robsinson et al. 2005 1. These UV observa- 
tions provided compelling evidence for two distinct classes of 
flares during this eruption, each characterized by significantly 
different ratios of their measured near UV (NUV) and far UV 
(FUV) fluxes. In this Paper we re-visit those data (which were 
analyzed under the assumption of a blackbody energy distri- 
bution) and in addition, we report on three further flare events 
on the stars CR Dra, AF Psc and SDSS J084425.9+513830.5 
(hereafter referred to as SDSS J084425) that were serendipi- 
tously recorded with a high time resolution of ~ 0.01 sec by the 
GALEX satellite during its all-sky survey observations during 
the period 2004-2005. 

2. Observations and Data Reduction 

All of the flare events reported in Table 1 were serendipitously 
recorded within the 1.2° field of view of the GALEX tele- 
scope in the far UV (FUV: 1350 - 1750A) and/or near UV 
(NUV: 1750 - 2800A) photometric bands during imaging sur- 
vey observations by the satellite (see Martin et al. 2005 for a 
detailed description of the GALEX instrument and its opera- 
tional modes). These observations are listed as either Medium 
Imaging Survey (MIS) or Deep Imaging Survey (DIS) modes 
in Table 1, together with their total on-orbit exposure times and 
their data-file identifier as listed in the Multi-Mission Archive 
at the Space Telescope Science Institute (MAST). The detec- 
tion of the flare on the star AF Psc was recorded during a 
pointed observation of a star-field selected for an observation 
of a galaxy under the auspices of NASA's GALEX Guest 
Investigator (GI) program. 

The FUV and/or NUV data for each flare were recorded 
as time-tagged events (with a time resolution of ~ 0.01 sec) 
using the GALEX photon counting detectors (Jelinsky et al. 
2003 1. These observations have been processed (as a non- 
standard data-product) using Version 4.0 of the GALEX Data 
Analysis Pipeline operated at the Caltech Science Operations 



Center (Pasadena, CA), which has been described in detail by 
Morrissey et al. ( 2005 1. The final data product is a flat-field cor- 
rected photometric time sequence of photons mapped in Right 
Ascension and Declination to the sky. Photon events associated 
with a stellar flare were summed within a 10 arcsec aperture 
centered at the expected position of the associated dMe star, 
while photons associated with the background emission were 
summed within an annulus extending from ~ 15 to 20 arcsec 
from the central position of the source. We assumed that any 
(faint) astronomical sources contained within the background 
annulus remained constant during the bright flare events. 

The photon count rates associated with each of the four 
flare events have been converted into UV fluxes of erg cm -2 
s _1 A -1 using the appropriate GALEX FUV and NUV-band 
calibrations listed in Morrissey et al. (2005). Using these flux 
corrected data we show the FUV and NUV light-curves as a 
function of elapsed time for the four flare events in Figure 1 
(we note that only NUV data was recorded for the flares de- 
tected on CR Dra and SDSS J084425 due to the FUV detector 
being turned off for satellite operational reasons). In Figure 2 
we also show the ratio of the measured FUV to NUV fluxes as a 
function of elapsed time for both AF Psc and GJ 3685 A, which 
(as described later) can be used as a probe of the temperature 
variation throughout the duration of these flares. 

A general comparison of the light-curves shown in Figure 
1 reveals that in all cases there was a very steep and rapid rise 
in flux (in both the NUV and FUV bands) over a time-scale of 
~ 70 seconds. Following this flare peak intensity all the flares 
exhibit an exponential-like decay that lasts for several hundred 
seconds. In the case of GJ 3685A there are multiple, smaller in- 
tensity flares during this decay phase. The smaller FUV flare on 
GJ 3685A, which is the least energetic of the events we have 
detected, is somewhat anomalous in that it seems to be com- 
posed of several closely spaced (in time) small flares that grad- 
ually increase in intensity over a time-span of ~ 40 seconds. 
The peak FUV flux then dies away very quickly to the pre-flare 
basal intensity in less than 50 seconds. Thus the character of 
the small flare on GJ 3685A exhibits an almost inverse behav- 
ior to the larger FUV and NUV flares recorded on the other 
stars, which generally consist of a very rapid increase in flux 
followed by an extended exponential-like flux decrease that is 
often associated with small flares prior to returning to the pre- 
flare background intensity. We also note that the NUV flare on 
CR Dra, which exhibits several small brightenings during its 
decay phase, does not exactly temporally track the FUV be- 
havior, suggesting an origin in a slightly different temperature 
region. 
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Fig. 1. The calibrated flare light-curves for AF Psc, CR Dra, GJ 3685A and SDSS J084425 for both the FUV and NUV bands. 



3. Data Analysis 

3. 1. Micro -flare activity 



Lu & Hamilton (1991) suggested that solar flares consists of 
an avalanche of many very small reconnection events termed 
micro-flares (or nano-flares). These small-scale phenomena are 
also thought to be of fundamental importance in the heating of 
the solar (and stellar) coronae (Parker 1988 Doyle & Butler 
I1985> . The statistical analysis of rapid time-sequence UV ob- 
servations of flare activity on two dMe stars by Robinson et 
al. i 19951 1199 9 1 has shown that over a time-scale of a few 
hours most of the observed increase in integrated flux arises 
in a slowly varying stellar background whose short-term vari- 



ability seems consistent with the presence of many overlapping 
micro-flare events. In the case of the very large flare observed 
by GALEX on GJ 3685A, the stellar UV flux increased by a 
factor of 5 in the pre-flare period prior to the main flare event 
(Robinson et al. 2005 ). However, a statistical analysis of the 
pre-flare activity of this energetic eruption failed to reveal com- 
pelling evidence for an increased level of micro-flaring events. 

Due to the relatively short duration (~ 1500 sec) of each 
flare observation by GALEX, the present time-tagged photon 
data are not ideally suited for a rigorous statistical investiga- 
tion of micro-flare activity, which is best performed over far 
longer periods (Robinson et al. 1995 1999 1. However, for the 
700 seconds of 'quiescent' pre-flare NUV data for CR Dra and 
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SDSS J084425 we have performed a simple test of the statis- 
tical significance of the variability in the recorded photon flux 
(Robinson et al. [1995 1. Briefly, this method involves search- 
ing the pre-fiare photon data for time periods during which 
the count rate deviates significantly above that expected from 
Poisson statistics. Large micro-flares reveal themselves if the 
data is binned into small time intervals (< 1 sec), whereas weak 
micro-flare events require binning over far longer periods (~ 
10 sec). A micro-flare is deemed a real event if it is detected 
with a probability < 10~ 2 of it occurring by chance (i.e. the 
microflare event is > 3er above the statistical mean of the data). 

In Figure 3 we show 700 seconds of the pre-fiare photon 
data for CR Dra and SDSS J0084425 binned into 0.2, 1.0 and 
10.0 second intervals together with the statistical mean levels 
and their respective ±3cr limits. Unfortunately none of these 
plots reveal the presence of a significant flux contribution from 
any major micro-flare activity. However, due to the limited pe- 
riod of pre-fiare observation our present findings do not pre- 
clude the presence of micro-flaring immediately prior to the 
beginning of our observations. 



DEM, we have derived the emission line spectrum over the 
wavelength interval 1200 to 3000 A. Multiplying by the FUV 
and NUV response functions then gives us the emission lines 
contained within the two filters. We note that the continuum 
contribution from two-photon, bremsstrahlung and free-bound 
processes are all included in these calculations. The FUV and 
NUV filter contributions were calculated for a range of elec- 
tron pressures from 10 15 to 10 19 cm~ 3 K. For pressures around 
10 16 - 10 17 cm~ 3 K, the resulting FUV:NUV filter ratios were 
around 0.2, substantially less than the observed quiescent value 
of -0.6 for AF Psc and -1 for GJ 3685 A. 



4. Results 

4. 1. The UV Flare Spectrum 

The GALEX far ultra-violet band response extends from 
-1350A to 1750A while the near ultra-violet band extends 
from -1750A to 2800A. As already seen from IUE, FUSE 
and HST spectra of active dMe stars, these spectral windows 
contain many emission lines. Based on data from the above sci- 
ence missions, the differential emission measure (DEM) distri- 
bution of dMe stars has been shown to be similar to the Sun. 
However, the lack of simultaneous observations over a suffi- 
ciently wide wavelength range has resulted in DEM distribu- 
tions derived over limited temperature intervals (del Zanna et 
al. I2002I Osten et al. I2005l >. In the present work we use the 
volume emission measure of Osten et al. (2005 1 as a starting 
point for the active dMe star EV Lac, and then convert it to 
the differential emission measure distribution (for a discussion 
of the different definitions see del Zanna et al. 12002b . Using 
the CHIANTI package (Landi et al. 120051 and the EV Lac 



Fig. 4. Sample DEM's for a dMe star in quiescence (Q) and 
flaring state (F) (see text) 



As mentioned above, several temperature intervals within 
the DEM distribution are poorly constrained, e.g. the upper 
chromosphere/lower transition region. Reducing the EV Lac 
DEM by a factor of five in the interval log T e = 4.0 to 4.6 re- 
sults in a FUVNUV ratio of 0.5 to 1.0. Furthermore, reducing 
the coronal DEM contribution around 2 x 10 6 K by a factor of 

gives FUVNUV ratios close to unity even at low electron 
pressures. In Table 2 we tabulate the FUVNUV filter ratios 
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Fig. 3. Pre-flare quiescent photon data for CR Dra (left) and SDSS J0084425 (right) binned in 0.2, 1 .0 and 10.0 second intervals, 
shown together with the 3-cr statistical error limits (dashed lines). No evidence is found for micro-flare activity. 



for such an amended quiescent DEM distribution as shown in 
Figure 4. 

For the flaring DEM distribution we increased the emis- 
sion measure over all temperatures ranges, in particular 
in the region around 100,000 K where the important C iv 
line is formed. In the third column of Table 2 we list the 
FUV:NUV filter ratios, here we see values above 3 for 
the high electron pressure. This shows that in addition to 
changing the DEM distribution, the electron pressure plays 
an important role. In earlier work on this data, Robinson 
et al. (20051 used blackbody radiation emission in an at- 
tempt to reproduce the observed measured FUV:NUV ra- 
tios, but found that an additional source of emission was 
required which he suggested was line emission. With a suit- 
able DEM distribution we derive FUV:NUV ratios ranging 

from 0.59 to 0.93 for electron pressures of 10 15 10 16 

cm 3 K in excellent agreement with the observed pre-flare 
values for AF Psc and GJ 3685A. Using the flare DEM and 
a higher electron pressure gives FUV:NUV ratios around 3, 



in good agreement with the AF Psc flare but a factor of two 
less than that observed in the GJ 3685A flare. Combining 
the above emission line contribution with a possible "white 
light" contribution resulting from a higher temperature 
than the star's photosphere could lead to larger FUV:NUV 
for the more energetic flares. Whether this "white light" 
contribution is important in the far UV remains uncertain. 

Regarding other emission contribution possibilities, 
previous studies of UV and X-ray emission in solar flares 
have shown an excellent temporal correlation (Cheng et 
al. [|981 1, indicating that they might both be a result of 
the same process of bombardment of the atmosphere by 
high energy particle beams. However, more recent data 
(Alexander & Coyner 120061 caution against such a pre- 
conception. These authors found that the UV and hard X- 
ray sources were spatially separated, and although they 
showed a high temporal correlation there was evidence for 
a magnetic connection between the two regions. Whether 
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this could enhance the UV continuum emission during dMe 
flares is unclear. 

A more likely continuum enhancement mechanism is 
silicon reconnection. Phillips et al. 1 1992 1 discussed the con- 
tinuum observed in the far-ultraviolet region from IUE 
spectra of dMe stars in a flaring state compared with that 
observed in solar flares. They suggested that the excita- 
tion of this continuum was due to the ionization of neu- 
tral silicon atoms near the temperature minimum region 
irradiated by ultraviolet line radiation emitted by the up- 
per chromosphere or transition region. This also appears 
to be the case for two flares on the RS CVn star (II Peg), 
and one on the spotted active binary BY Dra. The evidence 
is an observed proportionality of the continuum intensity 
with the intensities of strong transition region lines such 
as C iv 1548.51A and C n 1335.36A. This process is not 
presently included in CHIANTI and will affect the FUV fil- 
ter at wavelengths less than 1682A. 

Using our presently proposed flare DEM at an electron 
pressure of 10 18 cm~ 3 K, the dominant emitter is C IV, con- 
tributing 35% of the FUV filter with Si IV contributing another 
8%. The contribution from weaker lines and the continuum are 
not insignificant, e.g. between 1560 A and 1750 A, there is 
~27% contribution. In the NUV filter, the dominant line emit- 
ters (originating mostly in the upper chromosphere/lower tran- 
sition region) are Mg II, Fe II, Al III, C III, etc. For example, 
Mg II has — 10% contribution, the region from 2320 to 2530 A 
is -17% (mostly Fe II), while the lines around 2290 ± 20 A 
supplies —14%. The remaining NUV flux contribution arises 
from the various continuum processes. 

Several coronal lines are also present in both filters, e.g. 
Fe xxi 1354, Fe xi 1467, vn 1623, Fe xxi 2298, Fe xx 2666 
and Ca XVI 273 1 . The accumulated contribution of these coro- 
nal lines in the FUV filter is around 2%, rising to over 10% in 
the NUV filter. 

These calculations therefore suggest that changing the 
DEM distribution can change the FUVNUV flux ratio by a fac- 
tor of 2-3. However, a major influence on the ratio concerns the 
increase in electron pressure during the flare. Obviously, with a 
different DEM these derived emission fractions will change. In 
Figure 5 we show the range of possible spectral lines present in 
the GALEX FUV and NUV bands during flares on dMe stars 
based on the DEM (F) distribution from Figure 4. 

Very few M dwarfs have directly measured electron densi- 
ties in the transition region or corona. Testa et al. (2004), Ness 
et al. (2004) and Osten et al. (2005) give electron densities 
measurements for a few dMe stars including EV Lac. Using 
the O VII lines formed around 2.5 x 10 6 K, Testa et al. (2004) 
give N e — 5.6 x 10 10 cm~ 3 , rising to 5.6 x 10 12 cm~ 3 from 
Mg XI formed around 8 x 10 6 K. Values given by Osten et al. 
are slightly larger, while those given by Ness et al. are similar at 
2.5 x 10 6 K rising to just over 10 11 cm~ 3 at 4x 10 6 Kas derived 
from Ne IX. At lower temperatures, e.g. 2.5 x 10 5 K, Osten et 
al. derived an electron density of — 10 11 cm~ 3 from O V. This 
suggests electron pressures ranging from — 3xl0 16 cm -3 K 
at 2.5 x 10 5 K to - 4xl0 19 cm" 3 K at 8 x 10 6 K. It should 
however be noted that the derived density at these higher tem- 
peratures are a little uncertain due to line blending problems. 



Table 2. A summary of the GALEX FUV:NUV band flux ra- 
tios for the two different DEM distribution and five electron 
pressures as discussed in the text. 

electron pressure (cm" 3 K) FUVNUV FUVNUV 

DEM (Q) DEM (F) 



10 ib 0.59 1.00 

10 16 0.93 1.66 

10 17 1.27 2.35 

10 18 1.46 2.54 

10 19 1.84 3.04 



As an alternative to the contact pressure assumption, we re- 
run the models for a constant electron density of 10 11 cm -3 
and 10 12 cm -3 . For the DEM (F) distribution, this gives an 
FUV:NUV flux ratio of around 2.5. 

4.2. Oscillations 

Magneto-hydrodynamic oscillations associated with coronal 
loops in active regions on the Sun have been widely reported 
(Aschwanden et al. I1999t . but only recently have such oscil- 
lations been observed during an X-ray flare on an active M- 
dwarf star (Mitra-Kraev et al. 2005b). In this latter case ev- 
idence was found for a damped oscillation with a period of 
— 750 s, from which a coronal loop length of 2.5 x 10 10 cm 
was derived. Based on this recent success in finding oscilla- 
tions in M-dwarf flare spectra, we have similarly searched for 
periodicities in the observed UV photon count-rates from the 
four flare events observed by GALEX using two different an- 
alytic tools: the wavelet and the Fourier transform methods. 
The wavelet transform is used here since it provides informa- 
tion on the location in time of where the oscillations occur, their 
temporal spread, and also the period of the oscillations. In this, 
the use of wavelets are an advantage over the more traditional 
Fourier power spectra method, which only gives information 
on the periods present within a time series, but no informa- 
tion on where in time that period occurs within the time series 
of data. Details on the wavelet analysis, which provides infor- 
mation on the temporal dependence of a signal, are described 
in Torrence & Compo ( 1998 ) 1 . For the convolution of the time 
series in the wavelet transform we have chosen the Morlet func- 
tion as defined in Torrence & Compo ( 1998 ). We note that the 
wavelet transform "inspects" the time series at a number of 
temporal scales (periods). In effect, the wavelet transform is 
a bandpass filter (Torrence & Compo 1998 ), and so, by choos- 
ing an appropriate scale it is possible to filter the original time 
series. This is possible since the wavelet transform is a band- 
pass filter with a known response function and so it is possible 
to "reconstruct" the original time series using an inverse filter. 
In our case, we chose only those scales less than 60 sec and 
then reconstructed the time series containing only periods less 
than this. The results of applying this method can be seen in 
Figures 6 to 1 1 . To determine whether or not any oscillations 
that were found were real, we implemented the Linnell Nemec 
& Nemec ( 1985 ) randomization method that estimates the sig- 

1 see http://atoc.colorado.edu/research/wavelets/ 
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Fig. 5. Sample dMe star flare emission spectrum with line identifications after folding the spectrum derived from the (F) DEM in 
Figure 4 with an electron pressure of 10 18 cm -3 K through the GALEX FUV and NUV band response functions (see text). 
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Fig. 6. The left two plots show the flare light-curve plus the de-trended NUV data for AF Psc, while the central two panels sho 
the detected oscillation periods less than 60 s plus their associated probability. The right panel shows the result of the Fourier 
analysis. 



nificance level of the main peak in the wavelet spectrum. The 
use of the randomization technique was performed according 
to O'Shea et al. J200U . Only those periods with a significance 
level above 95% were considered real and used in this work. 
The use of a significance level of 95% underlines the reliability 
of the wavelet technique. To back up the use of the wavelets 
in detecting the oscillation periods, we also show estimates of 
periods obtained from power spectra measured using smoothed 
Fourier transforms. We used the Fourier transform method of 
Jenkins & Watts 
in Doyle et al. 



1968 1 in a similar manner to that outlined 
Significance levels of 99.9% were used 



19991 

for the power spectra displayed in Figures 6 to 11, and shown 
as the dashed lines. The fact that the wavelet and Fourier re- 
sults show the same periods is further evidence of the relia- 
bility of using wavelets to measure the periods of oscillation. 
We note that Figures 6-1 1 show the flare light-curve, the de- 



trended data eliminating all periods greater than 60 seconds, the 
wavelet probability level, the detected periods from the wavelet 
and Fourier results for all 4 stellar flares. The wavelet probabil- 
ity refers to the maximum power at each time period in the 
wavelet plot, and is indicated by the over-plotted white line. 

Finally, we note that using the "post-mortem" technique of 
Schwazenberg-Czerny ( 1991 1 and the power spectra in Figures 
6 to 1 1 we can estimate an error in the periods measured from 
the power spectra of approx. 3-5 sec, which in each case is 
much less than the typical periods of 30-40 sec being mea- 
sured. 

Evidence for short period oscillations (less than 60s) is 
found for all of the 4 flaring sources shown in Figure 1. In 
the NUV and FUV bands of AF Psc strong oscillations with 
a r^30s period are seen in the wavelet plots. This is confirmed 
by the results from the accompanying power spectrum where 
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Fig. 7. The left two plots show the flare light-curve plus the de-trended FUV data for AF Psc, while the central two panels sh< 
the detected oscillation periods less than 60 s plus their associated probability. The right panel shows the result of the Fourier 
analysis. 
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Fig. 8. The left two plots show the flare light-curve plus the de-trended NU V data for GJ 3685 A, while the central two panel™ 
show the detected oscillation periods less than 60 s plus their associated probability. The right panel shows the result of the 
Fourier analysis. 



broad peaks at ^30 sec are seen. In these plots of power spec- 
tra the large amount of low period (high frequency) noise can 
be easily seen. Similarly, both the NUV and FUV bands of 
GJ 3685A show evidence for short oscillation periods of ^40 
sec. The main peaks in the power spectra (below 60 sec) for 
these two bands are at the same period of ^40 sec. However 
the FUV band shows a more complicated structure with evi- 
dence of many more peaks present down to lower periods than 
that seen in the NUV band. For the NUV data of GJ 3685A, 
the wavelet analysis suggest periods of around 30 sec during 



the initial stages of the second flare, increasing to 40 sec dur- 
ing the decay phase. This is also seen in the Fourier results. 
Similarly, the power spectrum CR Dra NUV shows strong evi- 
dence for a period around 40 sec. In the wavelet spectrum plot 
for this source the oscillation is shown as a single peak at a pe- 
riod of ^47 sec, lasting for only the time period between 1 3 and 
16 minutes. The failure of the wavelet spectrum to resolve two 
peaks at ~40 and 50 seconds is related to its relatively poorer 
temporal resolution compared to the Fourier power spectrum, 
i.e., the restricted number of scales that can be chosen. This is 



10 



Welsh et al.: GALEX dMe flare events 



10000 
8000 
6000 
4000 
2000 



GJ 3685A 
FUV filter 



Level of Probability 

f- 




Time (minutes) 



Trend subtracted 



L J 



3x10" 



1x10* 




10 20 30 40 50 60 70 
Time (sec) 



16 

Time (minutes) . 
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Fig. 10. The left two plots show the flare light-curve plus the de-trended NUV data for CR Dra, while the central two panels show 
the detected oscillation periods less than 60 s plus their associated probability. The right panel shows the result of the Fourieft 
analysis. 



an inherent drawback of the wavelet method that is compen- 
sated by the information the wavelet gives us on where an os- 
cillation is occurring in time during the observing period. For 
the source SDSS J084425, in the NUV we again find a period 
of ^-40 sec, measured from both the Fourier and wavelet meth- 
ods. 



flare maximum, the above authors equated these two quantiles, 
giving the loop length as 

1500 



L 



(1 _ 3.1.58)4/7 •d 



4/7 3/7 



(1) 



Several authors have used radiative and conductive cooling 
times (e.g., Hawley et al. 1995 1 to estimate flare loop lengths. 
During the rise phase, strong evaporative heating is dominant 
while during the decay phase radiative cooling is dominant. At 



where is the flare decay time, r,. is the flare rise time, T the 
apex flare temperature, and Xd — \J{cdl c max ), with c max the 
peak flare count rate and Cd the count rate at the end of the 
flare. For the present flares, the typical rise time is ~ 70 sec, 
the decay time ~ 300 sec and Cd/c max ~ 0.1. Unfortunately, 
we do not know the plasma temperature, but assuming a range 
of temperatures from 5 x 10 6 K to 20 x 10 6 K, implies loop 
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Fig. 11. The left two plots show the flare light-curve plus the de-trended NUV data for SDSS J084425, while the central two 
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lengths of 0.5 — 1 x 10 9 cm. This temperature range is in broad 
agreement with the results from observations, e.g. in the mod- 
elling of X-ray data, Sciortino et al. ( 1 19991 1 derived a maxi- 
mum temperature of 20 x 10 6 K for EV Lac and a much higher 
temperature of 300 x 10 6 K for AD Leo. Several authors have 
also modelled the non-flare X-ray spectral data, e.g. Osten et 
al.( 2005 1 and Sciortino et al. dl999> used two and even a three- 
component model to fit their data. The implied loop lengths 
from their work was 0. 1R* for loop temperatures around a few 
millions degrees, with even smaller loops at the higher temper- 
atures. 

Based on the spectral analysis outlined in the previous sec- 
tion, the implied electron pressure during a flare may be close 
to 10 18 — 10 19 cm -3 K, implying large electron densities, con- £ 



sistent with measurements. Emslie (1981) has suggested that 
quasi-periodic increases in the hard X-ray flux from solar flares 
may be the result of the local gas pressure attaining very large 
values, significantly larger than the values suggested in the pre- 
vious section. However, the idea that transient increases in the 
electron density result in quasi-periodic increases in the line 
flux is never-the-less valid. The quasi-periodic increases in the 
line flux could be interpreted in terms of fundamental mode 
quasi-standing waves, these being the easiest to excite. Such 
interpretation has already been given by Mitra-Kraev et al. 
( 2005b ) to flare oscillations on an active M-dwarf star. The in- 
terpretation in terms of propagating waves or higher harmonics 
is less plausible given the spatial resolution of the observations 
and the cancellation effect along the loop, i.e., any density in- 
crease on one side of the loop would be accompanied by a cor- 
responding decrease one the other side. The theory of coronal 
loop oscillations was developed by Roberts et al. {1984). The 
periods of the oscillations are given by the expression 



where c corresponds to either the fast kink mode, the fast 
sausage mode, or the slow sausage mode. The fast kink mode 
is almost incompressible and can therefore be discarded from 
further discussion. According to Roberts et al. ( 1984 >, the 
wavenumbers required for the existence of the fast sausage 
wave are high under coronal conditions due to the cutoff prop- 
erty. In the case of a standing mode, the wavenumber is deter- 
mined by the loop dimensions and is below the cutoff value. 
The remaining possibility is the slow sausage mode. The rela- 
tionship between the period P, loop length L and temperature T 
for the fundamental mode slow standing wave is approximately 
given by 



P 



c T or L = 7600Pc Ts /T, 



(3) 



with cts being the ratio between the cusp speed ct and the 
sound speed c s . Comparing the above relationship with for- 
mula Q we find a ratio ct s = 0.73 for a period of 40 s. In com- 
parison, Mitra-Kraev et al. J2005bt found a ratio ct s = 0.82. 
This means that the plasma inside the loop (i.e., the ratio 
between the thermal pressure and the magnetic pressure) is 
slightly higher in the present case. 

Slow sausage (acoustic) fundamental mode standing waves 
have been recently detected on the Sun (Wang et al. 2003 1. 
These waves are caused by the injection of hot chromospheric 



P = 



2L 



(2) 



plasma at the footpoints of the loops. Taroyan et al. (|2005l 
showed that such oscillations could be readily excited if the 
duration of the plasma injection at the loop footpoint matches 
the natural fundamental mode period of the loop. The closeness 
of these two quantities would result in quasi-standing acoustic 
oscillations. In the present case, the injection of plasma into the 
loop could be associated with the flare and should last approx- 
imately 40 s. Finally, on a precautionary note, we acknowl- 
edge that the sequential heating of evenly spaced loops at 
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a constant rate (i.e. in an arcade) could also reproduce the 
oscillation signature that we have found in these data. We 
refer the reader to the excellent review of coronal oscilla- 
tions by Roberts 1 2000 1 for further details of this, and other 
possible oscillation mechanisms. 

5. Conclusion 

We have presented near and far ultraviolet light-curves for stel- 
lar flares observed with the NASA GALEX satellite on the 
four nearby dMe-type stars GJ 3685A, CR Dra, AF Psc and 
SDSS J084425. These data were recorded as time-tagged UV 
photon events with a temporal resolution of < 0.01 seconds, 
thus enabling a detailed time-analysis of both their quiescent 
and flaring states. During 700 seconds of quiescence, prior 
to the flare outbursts on both CR Dra and SDSS J084425, 
no statistical evidence was found for micro-flare outbursts in 
the NUV with data binned in 0.2, 1.0 and 10.0 second time- 
intervals. 

We have used a modified differential emission model 
(DEM) curve of the flare star EV Lac in conjunction with the 
CHIANTI atomic database to produce theoretical DEM curves 
for a typical dMe star in both quiescence and a flaring state. 
From the quiescent DEM distribution we were able to repro- 
duce FUV:NUV flux ratios greater than unity, consistent with 
the measurements during the flares. A critical input parameter 
in determining the FUV:NUV ratio is the value of plasma elec- 
tron density during the flare event. It was also found that the 
major single emitter in the FUV filter is C IV, although weaker 
emission lines and the continuum provide ^50% of the flux. 
In the NUV filter, a significant contribution comes from upper 
chromosphere/lower transition region emission lines, in addi- 
tion to continuum emission. Several coronal lines are present in 
both filters, with ~10% contribution in the NUV filter, whose 
main contributer is from continuum processes. 

We have also searched for periodicities in the observed 
UV count-rate data recorded during the four flare events to 
search for coronal loop oscillations. Using both the wavelet and 
Fourier transform analysis methods on these flare data, we have 
found compelling statistically significant evidence for oscilla- 
tions in both of the NUV and FUV data for flares on AF Psc 
and GJ 3685A, and in the NUV data for the flares on CR Dra 
and SDSS J084425 with periodicities in 30-40 second range. 
Assuming typical values of electron pressure during a flare, 
these periodicities may be due to acoustic waves in coronal 
loops of length ~ 10 9 cm for a plasma with a temperature range 
of 5 — 20 x 10 6 K. These values suggest a loop length of less 
than 1/10" 1 of the M-dwarf stellar radii. Finally we note that 
although oscillations have been previously observed on the Sun 
and in X-ray observations of the M-type dwarf AT Mic, we be- 
lieve that this is the first detection of non-solar coronal loop 
flare oscillations at ultraviolet wavelengths. 
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